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ABSTRACT 

After at least 6 years of quiescence, Anomalous X-ray Pulsar (AXP) 4U 0142+61 entered an active 
phase in 2006 March that lasted several months and included six X-ray bursts as well as many 
changes in the persistent X-ray emission. The bursts, the first seen from this AXP in >11 years 
of Rossi X-ray Timing Explorer monitoring, all occurred in the interval between 2006 April 6 and 
2007 February 7. The burst durations ranged from 0.4— 1.8xl0 3 s. The first five burst spectra are 
well modeled by blackbodies, with temperatures kT ^ 2 — 9 keV. However, the sixth burst had a 
complicated spectrum that is well characterized by a blackbody plus two emission features whose 
amplitude varied throughout the burst. The most prominent feature was at 14.0 keV. Upon entry 
into the active phase the pulsar showed a significant change in pulse morphology and a likely timing 
glitch. The glitch had a total frequency jump of (1.9±0.4)xl0 -7 Hz, which recovered with a decay 
time of 17+2 days by more than the initial jump, implying a net spin-down of the pulsar. Within 
the framework of the magnetar model, the net spin-down of the star could be explained by regions of 
the superfluid that rotate. slower than the rest. The bursts, flux enhancements, and pulse morphology 
changes can be explained as arising from crustal deformations due to stresses imposed by the highly 
twisted internal magnetic field. However, unlike other AXP outbursts, we cannot account for a major 
twist being implanted in the magnetosphere. 

Subject headings: — stars: neutron — X-rays: stars — X-rays: bursts pulsars: individual 
(4U 0142+61) 


1. INTRODUCTION 

It is now generally accepted that the class of objects 
referred to as “Anomalous X-ray Pulsars” (AXPs) are 
magnetars - young isolated neutron stars powered by 
the evolution of their high magnetic fields (Thompson &: 
Duncan 1995; Thompson & Duncan 1996). AXPs are 
X-ray pulsars with periods in the range 2-12 s, and X- 
ray luminosities (~ 10 33 — 10 35 erg s~^) that cannot be 
accounted for by their available spin-down energy. The 
magnetar model was first proposed to explain the dra- 
matic behavior exhibited by an apparently different ob- 
ject class - the Soft Gamma Repeaters (SGRs). SGRs 
show persistent properties similar to AXPs, but they 
were first discovered by their enormous bursts of soft 
gamma rays (> 10 44 erg) and their much more frequent, 
shorter, and thus less energetic bursts of hard X-rays. 
To date, SGR-like X-ray bursts have been observed from 
six AXPs, thus solidifying the connection between the 
two source classes (Gavriil et al. 2002; Kaspi et al. 2003; 
Woods et al. 2005; Krimm et al. 2006; Mereghetti et a. 
2009; Kumar & Safi-Harb 2011). For reviews of magnetar 
candidates and AXPs see Woods & Thompson (2006), 
Kaspi (2007), and Mereghetti (2008). 

Thus far, only the magnetar model can explain the 
bursts observed from SGRs and AXPs (Thompson & 
Duncan 1995). The internal magnetic field exerts stresses 
on the crust which can lead to large scale rearrangements 
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of the external field, which we observe as giant flares. If 
the stress is more localized, then it can fracture the crust 
and displace the footpoints of the external magnetic field 
which results in short X-ray bursts. The highly twisted 
internal magnetic field also slowly twists up the exter- 
nal field; magnetospheres of magnetars may therefore be 
globally twisted (Thompson et al. 2002). Reconnection 
in this twisted magnetosphere has also been proposed as 
an additional mechanism for the short bursts (Lyutikov 
2002). 

In addition to bursts, AXPs and SGRs exhibit pulsed 
and persistent flux variations on multiple different 
timescales. An hours-long increase in the pulsed flux 
has been seen to follow a burst in AXP IE 1048.1-5937 
(Gavriil et al. 2006). On longer timescales, AXPs can 
exhibit abrupt increases in flux which decay on several- 
week timescales. These occur in conjunction with bursts 
and have been suggested as being due to thermal ra- 
diation from the stellar surface after the deposition of 
heat from bursts. Such flux enhancements have been 
observed in SGRs (see Woods et al. 2001, for exam- 
ple). The flux enhancement of AXP IE 2259+586 dur- 
ing its 2002 outburst can also interpreted as burst af- 
terglow (Woods et al. 2004), however, a magnetospheric 
interpretation is also plausible (Zhu et al. 2008). AXP 
IE 1048.1-5937 exhibited three unusual flux ‘flares.’ In 
the first two, the pulsed flux rose on week-long timescales 
and subsequently decayed back on time scales of months 
(Gavriil & Kaspi 2004; Tam et al. 2008). These vari- 
ations have been tentatively attributed to twists im- 
planted in the external magnetosphere from stresses on 
the crust imposed by the internal magnetic field. AXPs 
XTE J1810-197, IE 1547.0-5408, and the AXP candi- 
date AX J1845-0258 have also exhibited large flux vari- 
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in IE 2259+586 or the slow-rise type as in IE 1048.1- 
5937. CXO J164710. 2-455216 showed a clear abrupt 
rise (Muno et al. 2007). 1RXS J1708-4009 has been ar- 
gued to have flux variations associated with timing events 
(Israel et al. 2007b). Finally, AXP 4U 0142+61 has ex- 
hibited the longest timescale flux variations thus far, in 
which the pulsed flux increased by 29+8% over a period 
of 2.6 years (Dib et al. 2007; Gonzalez et al. 2010). 

4U 0142+61 is an 8.7-s AXP. It has a period derive 
tive of P = 0.2 x 10“ n , implying a surface dipole mag- 
netic field of 1.3 x 10 14 G. 4U 0142+61 was monitored 
by RXTE in 1997 and from 2000 to 2007. Gavriil & 
Kaspi (2002) showed that 4U 0142+61 generally rotates 
with high stability. Morii, Kawai, & Shibazaki (2005) 
reported a possible timing glitch in 1999 on the basis 
of an Advanced Satellite for Cosmology and Astrophysics 
( ASCA ) observation in which the value of the frequency 
was marginally discrepant with the frequency as reported 
by Gavriil & Kaspi (2002). Dib et al. (2007) showed that 
the glitch may have occurred but is not required by the 
existing data. Dib et al. (2007) also reported on the evo- 
lution of the properties of 4U 0142+61 from 2000 March 
to 2006 April. In particular they reported stable timing, 
and an evolution of the pulse profile in 2—4 keV where 
the dip between the two peak was rising between 2000 
and 2006. They also reported a 29+8% increase in the 
pulsed flux between 2002 May and 2004 December. As of 
2006 March, in the published flux history of this source, 
there had been no reports of any X-ray activity such as 
bursts or flares, as described above for other AXPs. 

Here we report on the first detection of bursts from 
AXP 4U 0142+61, making this the seventh AXP for 
which this phenomenon has been observed. We also re- 
port that the source appears to have entered an active 
phase in 2006 March in which almost every aspect of 
the emission changed. Our observations are described in 
Section 2. Our burst, pulsed morphology, pulse phase, 
pulsed flux, and timing analysis are presented, respec- 
tively, in Sections 3.1, 3.2, 3.3, 3.4, and 3.5. In Section 4, 
we discuss the possible origins of this behavior and the 
implications for the magnetar model. 

2. OBSERVATIONS 

4U 0142+61 has been monitored with Rossi X-ray 
Timing Explorer ( RXTE) in 1997 and from 2000 to 2007. 
It has been monitored bi-monthly since 2005 March, with 
a typical observation length of ~5 ks. On 2006 March 
23, the source entered an active phase where many pulsed 
flux, spectral, and pulse profile changes were observed. 
We detected 6 bursts in three observations after the en- 
try into the active phase. After each burst detection, 
several RXTE target-of-opportunity (ToO) observations 
were made in addition to the regular monitoring. 

Here, we present a detailed analysis of the three obser- 
vations containing bursts (see Table 1). We also present 
an analysis of 91 other observations spanning the 2005 
March 21 to 2007 May 15 time period (MJD 53481 to 
MJD 54235, Observation IDs 91070-05-04-00 to 92006- 
05-32-00). 32 of these observations were prior to the 
entry into the active phase, and the remaining 59 ob- 
servations were after. This long-term analysis was per- 
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Fig. 1. — Epochs of the RXTE observations of 4U 0142-f61 
analyzed in this paper (see § 3.2). Each observation is represented 
by a vertical line. The dashed line marks the entry of the source 
into the active phase. An M x” marks the first observation where the 
pulse profile was significantly different from the long-term average, 
“bl”, “b2”, and “b3” mark the observations containing bursts. To 
the left of the dotted line, the letters refer to groups of observations 
having similar total integration times. To the right, the letters refer 
to groups of observations having similar pulse profiles. 

formed in order to track the evolving pulsed flux and 
timing properties of the source. 

Figure 1 shows a timeline of the 94 analyzed observe 
tions. Note that two segments of any observation that 
was split counted as separate observations if the seg- 
ments were given different observation IDs. The ranges 
of epochs with an increased density of observations con- 
tain the ToO observations. Prior to the active phase, 
the groups of observations referred to by capital letters 
in the Figure have similar total integration times. The 
groups of observations in the active phase have similar 
pulse profiles. These groups will be referred to in Sec- 
tions 3.2 and 3.5. 

3. ANALYSIS AND RESULTS 

All data presented here are from the Proportional 
Counter Array (PCA; Jahoda et al. 2006) aboard 
RXTE. The PCA is made up of five identical and in- 
dependent proportional counter units (PCUs). Each 
PCU is a Xenon/methane proportional counter with 
a propane veto layer. The data were collected in ei- 
ther GoodXenonwithPropane or GoodXenon mode which 
record photon arrival times with ~l-/us resolution and 
bins them with 256 spectral channels in the ~2-60 keV 
band. For all RXTE observations of 4U 0142+61, pho- 
ton arrival times at each epoch were adjusted to the solar 
system barycenter using the position obtained by Patel 
et al. (2003) from Chandra X-ray Observatory data. 

3.1. Burst Analysis 

For each monitoring observation of 4U 0142+61, using 
software that can handle the raw telemetry data, we gen- 
erated 31.25-ms lightcurves using all Xenon layers and 
events in the 2-20 keV band. These lightcurves were 
searched for bursts using our automated burst search al- 
gorithm introduced in Gavriil et al. (2002) and discussed 
further in Gavriil et al. (2004). In an observation on 2006 
April 6, we detected a significant burst, and four more 
bursts were detected in a single observation on 2006 June 
25. The sixth and most energetic burst was detected on 
2007 February 7. The bursts were significant in each ac- 
tive PCU. See Table 1 for the number of active PCUs in 
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servations Containing Bursts 


rvation ID 92006-05-03-00 92006-05-09-01 

Burst Date 6 April 2006 26 June 2006 

Burst MJD (UTC) 53831 53911 54 

Number of Active PCUs 3 a 3 

Number of Bursts 1 4 

Burst Temporal and Spectral Properties 


Burst Number 
t P c (UT) 

T r d, ‘(ir.s) 

Tf 0,1 (ms) 

Tt f> ' (s) 

f p g,i (counts/s/P CU) 
fp g,i,k (10" 10 erg/s/cm 2 ) 
ft h ’ i (counts/s/PCU) 
ft h ’ i,k (10" 10 erg/s/cm 2 ) 
Ftot j (counts/PCU) 

Ftot ,i,k (10” 10 erg/cm 2 ) 

T 90 1 (s) 

Phase m (cycles) 
kT n (keV) 

«°( m) 

v 2 /DoF p fDoFl 


B One PCU switched on and another one off part way through the observation, but the total number of active PCUs stayed constant 
throughout the observation. 

b All quoted errors represent 1 -u uncertainties. In the fits, the column density was held fixed at N’//=0.64xl0 22 cm" 2 (Durant & 
van Kerkwijk 2006b). 
c Time of burst peak. 
d Rise Time. 

e Short-term decay time of the tail of the burst. 
f Long-term decay time of the tail of the burst. 
g 2-60 kev peak flux. 

h 2-60 keV flux in the tail component of the burst. 

1 These parameters were obtained by fitting the model given by Eq. 1 to the burst. 

J 2-60 keV fluence. These values were obtained by using Eq. 2. 
k All quoted fluxes and fluences in CGS units are unabsorbed. 

1 T90 duration, defined as the time it took to collect 90% of the total burst fluence. These values were obtained by using Eq. 3. 
m The phases are relative to the template shown in Figure 7. 

n Blackbody temperature. For burst 6, the temperature is that of the blackbody component after accounting for the apparent spectral 
features (see §3.1.2). 

0 The blackbody radius was calculated assuming a distance to the source of 3.5 kpc (Durant & van Kerkwijk 2006a). 
p Reduced x 2 ^ the spectral fits. 
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each burst observation, as well as for the burst observa- 
tion epochs. 

To analyze these bursts we created event lists in FITS 4 
format using the standard FTODLS 5 . For consistency with 
previous analyses of SGR/AXP bursts we extracted 
events in the 2-60 keV band. The burst lightcurves are 
displayed in Figure 2 (LEFT). Our techniques for char- 
acterizing the temporal and spectral properties of bursts 
were discussed in detail in Gavriil et al. (2004) but we 
repeat them here as some required modification in order 
address the distinct properties of these bursts. 

Before measuring any burst parameters we deter- 
mined the instrumental background using the FT00L 
pcabackest. We extracted a background model 
lightcurve using the appropriate energy band and num- 
ber of PCUs. pcabackest determines the background 
count rate only on 16-s time intervals, so we interpolated 
these values by fitting a polynomial of order 6 to the en- 
tire observation, which yielded a good fit for each data 
set. 


3.1.1. Burst Temporal and Energetic Properties 

The burst peak time ( t p ) was determined using the 
method described in Gavriil et al. (2004); Gavriil et al. 
(2006). Usually, to measure the fluence for SGR and 
AXP bursts, we subtract the instrumental background 
for the lightcurve, integrate the light curve and fit the 
result to a step function with a linear term whose slope 
is the “local” background rate. The burst fluence in this 
case is the height of the step function. Although this 
technique worked well for the first burst, which was a 
short, isolated event, it was not appropriate for bursts 2, 

3, and 4 because they had overlapping tails, and bursts 

4, 5 and 6 had tails that extended beyond the end of the 
observation. Thus, we opted to determine the best-fit 
burst profile , /(t), consisting of exponential rises and 
decays, 




b(t) + B 


-\-b(t) + B 


4 http://fits.gsfc.nasa.gov 

6 http: / 'heasarc.gsfc.nasa.gov/docs/software/ftools/ 
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fixed at the values determined by the methods outlined 
above. We fit for the rise time (r r ), the peak flux (/ p ), 
the flux of the long-term tail component (/*), and for 
the short (t/) and long-term (r t ) decay times of the tail 
of the burst. Figure 2 (LEFT) shows our observed and 
model burst lightcurves. The binning was chosen such 
that the flux in the peak bin was equal to the flux de- 
termined from our maximum likelihood fits. We then 
integrated our events to obtain fluence lightcurves. Fig- 
ure 2 (RIGHT) shows our observed and model burst flu- 
ence time series. For both figures, the contribution from 
neighboring bursts has been subtracted. 

As is done for 7-ray bursts and SGR and AXP bursts, 
we characterized the burst duration by T90, the time 
when 90% of the total burst fluence has been collected. 
From Eq. 1 it follows that the total fluence, F to t, is given 
by 

Ftot = fpTr + ( fp — ft)Tf + frTt , ( 2 ) 

and the X90 duration is given by 

Too = r t log (10/tT t /Ftot) • (3) 

Although Eq. 2 is exact, Eq. 3 is an approximation that 
holds very well because T90 r r ,T/ for all bursts. This 
fact is made obvious by Figure 2 (RIGHT). All burst 
temporal parameters are presented in Table 1. 

3.1.2. Burst Spectral Properties 

Burst spectra were extracted using all the counts 
within their Tqo interval. We then filtered the spectra 
to include counts only from the Xenon layers and from 
all active PCUs other than PCU 0. PCU 0 was ignored 
because of the loss of its propane layer and due to its fre- 
quent breakdowns. Background intervals were extracted 
from long, hand-selected, featureless intervals prior to 
the bursts. Response matrices were created using the 
FT00L pcarsp. The burst spectra were grouped such 
that there were at least 10 counts per bin after back- 
ground subtraction. Burst spectra, background spectra, 
and response matrices were then read into the spectral 
fitting package XSPEC 6 vl2.3.1. The spectra were fit to 
photoelectrically absorbed blackbodies using the column 
density found by Durant & van Kerkwijk (2006b). Only 
bins in the 2-30 keV band were included in the fits. The 
blackbody model provided an adequate fit for bursts 1 
through 0; model parameters are presented in Table 1. 

Extracted over its very long T 90 interval (see Table 1), 
the spectrum of Burst 6 was well fit by a simple black- 
body, however, there were hints of spectral features. On 
much shorter timescales, the spectrum of burst 6 was 
not well modeled by any simple continuum model and 
the presence of spectral features is clear (see Fig. 3). We 
looked closely at the first 8.69 s (one pulse period) of the 
burst (see Fig. 3), and determined that a simple contin- 
uum model with less than two spectral features was not 
a statistically good representation of the data. A simple 
photoelectrically absorbed blackbody yielded a reduced 
X 2 of 2.1 for 39 degrees of freedom (DoF). Adding a sin- 
gle Gaussian emission lines improved the fit but not suf- 
ficiently, yielding reduced x 2 s of 1.3 for 36 DoF. Only 

6 http://xspec.gsfc.nasa.gov 
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at energies 8.6lo!2 > 14.21q J. Note, that 8.6 keV line was 
very narrow, we thus kept its width fixed at the detectors 
energy resolution, this was not necessary for the 14.2 keV 
line as it was relatively broad. Our results from fitting 
different models to the first 8.69 s of the burst spectrum 
are summarized in Figure 3, and the parameters returned 
from our best-fit model are listed in Table 2. 

We investigated whether different continuum models 
could describe the feature rich burst spectrum just as 
well as a blackbody. We tried a simple power-law model, 
as well as more complicated continuum models such as a 
Cutoff Power-law (PLCUT; White et al. 1983) and a Neg- 
ative/Positive Exponential (NPEX; Mihara 1995). The 
PLCUT and NPEX models are often used to fit the con- 
tinuum component of the line rich spectra of accreting 
X-ray pulsars (see Coburn et al. 2002). All three of our 
alternate continuum models also required two emission 
features at similar energies as the blackbody, however 
they did not provide improved fits. The best alternative 
to a blackbody was the PLCUT model with two Gaussian 
emission lines at 8.8 and 14.3 keV with reduced x 2 =1-1 
for 33 DoF. Notice how this model combination had the 
same reduced x 2 as the blackbody combination, but re- 
quired one more free parameter. 

We then determined whether the features could be de- 
scribed by absorption rather than emission. We tried 
three different absorption models: cyclotron absorption, 
Gaussian absorption (multiplicative model), and simply 
subtracting Gaussians. We tested all three of these ab- 
sorption models with all the continuum models described 
above. No permutation provided a better fit than a 
blackbody with two emission lines, a simple power-law 
with two Gaussian absorption lines at 10.8 keV and a 
broad (Gaussian width a =4.4 keV) feature at 2.5 keV 
being the best absorption model combination with re- 
duced x 2 — 1-2 for 33 DoF. 

To study the temporal evolution of these features we 
extracted 17.38-s (two pulse periods) long spectra in 
steps of 4.34 s (half a pulse period) from the peak of the 
burst, and repeated the spectral fitting procedure out- 
lined above. Fig. 4 displays the spectral evolution of the 
features. The middle panel is a surface plot of the change 
in x after subtracting the continuum model (blackbody). 
A vertical slice in the middle panel corresponds to a spec- 
trum of the burst extracted over a two-pulse-period long 
interval, and these spectra were extracted in time steps 
of half a spin-period. The features showed clear tempo- 
ral variability and were most prominent near the onset 
of the burst (see Fig. 4). From the surface plot in Fig. 4 
(middle panel), we see that the 14.2 keV is highly sig- 
nificant at the start of the burst, and remains detectable 
at the > 1 -a level for ~130 s. The significance of the 
8.6 keV feature is intermittent, but it remains detectable 
at the > 1 -<r level beyond 130 s from the burst onset. 

3.2. Pulse Profile Changes 

Many AXP outbursts are accompanied by significant 
pulse profile changes (e.g. Kaspi et al. 2003). To search 
for these in 4U 01424-61, for each observation, we gen- 
erated 64-bin pulse profiles using the method described 
in Dib et al. (2007). We then aligned the profiles with 
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Fig. 2. — LEFT - The gray filled curves represent the 2-60 keV burst lightcurves as observed by RXTE. The binning was chosen such 
that the flux in the peak is equal to the model peak flux, f p (see Eq. 1), as determined by maximum likelihood likelihood testing. This 
binning scheme allow us resolve the bursts at the finest permissible time resolution. In descending order, the time bin widths are: 11.0, 2.3, 
107.6, 18.3, 32.6, and 686.3 ms.The solid black curves are the best-fit exponential rise and decay models (Eq. 1). RIGHT - The gray filled 
curves represent the 2-60 keV integrated burst lightcurves as observed by RXTE. For each burst lightcurve we subtracted the instrumental 
background, the background due to other sources in the FOV, and the contribution of the tails of neighboring bursts before integrating. 
Notice that some bursts have very long tails. The solid black curve is the integral of the model in the left panel. The dashed curve is a 
step function whose height corresponds to the total flux as determined by Eq. 2. 


a high signal-to-noise template using a cross-correlation 
procedure in the Fourier domain. Then, for each group of 
observations in the active phase, we summed the aligned 
profiles, extracted the DC component from the summed 
profile, and scaled the resulting profile so that the value 
of the highest bin is unity and the lowest bin is zero. 

The resulting pulse profiles are shown in Figure 5 in 
chronological order. The panels in the left-most column 
are from Dib et al. (2007) and show pulse profiles in the 
6 years prior to the entry into the active phase. Notice 
the slow change in the height of the dip between the 
peaks. The panels in columns 2 and 3, marked with the 
corresponding letter in the top right corner, show pulse 
profiles for each of the data groups in the active phase 
that were marked with a letter in Figure 1. The 4 plots 
in bold correspond to the observations marked with “x” , 
“bl”, “b2” , and “b3” in Figure 1. 

The pulse profile evolution can be described as fol- 
lows: prior to the active phase, features in the double 
peaked pulse profile were evolving on a timescale of sev- 
eral years (see panels C6 to CIO corresponding to RXTE 
Cycles 6 to 10). Then, in the first observation of the ac- 
tive phase (panel x), the pulse profile became suddenly 
triple-peaked. It was also triple-peaked in the following 
observation (panel bl), in which a burst occurred. It re- 
mained multi-peaked for the following two groups of ob- 
servations (panels G and H). Then, an observation with 


multiple bursts occurred (panel b2). In that observation, 
also having a triple-peaked pulse profile, the middle peak 
was taller than the other two. Following this observation, 
the pulse profile seemed to be slowly recovering back to 
its double-peaked long-term shape (panels I, J, K, and 
L). Another burst observation interrupted this evolution 
seven months later (panel b3). In that observation, a 
large burst was detected. In the pulse profile of that ob- 
servation, the left-most peak was significantly taller than 
usual. The event that caused this change had apparently 
no effect on the following observation which occurred 2 
days later: the profile went back to being double-peaked 
(panels M and N). To summarize, the pulse profile be- 
came multi-peaked a.t the beginning of the active phase. 
Following the second observation with bursts, the profile 
started to recover to its double-peaked shape. The evo- 
lution was only temporarily interrupted for the duration 
of the third observation with bursts. The behavior of the 
pulse profile in the 2—4 and 4—10 keV bands was similar. 

Note that from Figure 5 alone we can compare the 
sizes of the peaks to each other, but we cannot track the 
evolution of the heights of each peak separately. In order 
to do that, we must scale the pulse profile of each group 
of observations by the average pulsed flux of that group. 
This analysis is presented in Section 3.4.3. 

We also performed an analysis of the Fourier compo- 
nents of the pulse profiles. The results are shown in Fig- 










Fig. 3. — Panel 1A: The spectrum of burst 6 extracted over a 
pulse period (8.69 s) starting 50 ms before the burst peak. The 
red solid lines represents the best-fit model- a photoelectrically 
absorbed blackbody (red curve). The fit had a reduced x 2 of 2.1 
for 39 degrees of freedom (DoF). Panel IB: The difference between 
the data and the model shown in Panel 1A, divided by the uncer- 
tainty. Panels 2A and 2B : Same, except here the best-fit model 
(red curve) consists of a blackbody (dashed blue curve) plus a 
Gaussian emission lines (dashed green curve), all photoelectrically 
absorbed. The fit had a reduced x 2 of 1-3 for 36 DoF. Panels 3A 
and 3B: Same, except here the best-fit model (red curve) consists 
of a blackbody (blue curve) plus a two Gaussian emission lines 
(dashed green curves), all photoelectrically absorbed. The fit had 
a reduced \ 2 of 1.1 for 34 DoF. 

ure 6. The variations in the power of the fundamental 
Fourier component are shown in panel 1, and that of the 
second harmonic in panel 2. Note how the amplitude of 
the fundamental varied monotonically prior to, but not 
during, the active phase. Also note how the power in the 
second harmonic was already back to its pre-active-phase 
level before the last burst occurred. 


Parameter 

Value a 

Column Density b , Nh (10 22 cm~ 2 ) 

0.64 

2-60 keV Flux 0 (10“ 10 erg s^cm*- 2 ) 

39-6l| l 

Reduced x 2 (Degrees of Freedom) 

1.1(34) 

Blackbody Component 


Temperature, kT (keV) 5. 8+^3 

Radius 0 (km) 0. 19^01 

2-60 keV Flux c (10~ 10 erg s _1 cm“ 2 ) 34.8^5 


Gaussian Emission Line 1 


Energy (keV) 

Width e , a (keV) 

2-60 keV Flux° (10“ 10 erg s^ 1 cm“ 2 ) 

8-6ioi 

o.6±g;l 

Gaussian Emission Line 2 

Energy (keV) 

Width, <j (keV) 

2-60 keV Flux° (10 -1 ° erg s^ 1 cm“ 2 ) 

14.2l°;| 

■1 o + 0.3 
i *‘ 5 -0.3 

^-0.8 


a All errors were extracted at the 1 <7 level. 

b In all spectral fits, the column density was held fixed at the value 
found by Durant & van Kerkwijk (2006b). 
c All quoted fluxes are unabsorbed. 

d The blackbody radius was calculated assuming a distance to the 
source of 3.5 kpc (Durant & van Kerkwijk 2006a). 
e The width of this Gaussian line was not resolved, thus a Delta 
function was used smeared by the instrument response. 

made of the 5 Fourier harmonics that best fit the his- 
togram. The bottom curve is made of the best-fit 5 har- 
monics after the removal of the 4 seconds centered on 
each burst. Note how the two curves in a given plot are 
similar, demonstrating that the additional peaks in the 
profiles are not due to the burst. In fact, for the brightest 
(and longest) burst, a fold of the last 1700 seconds of the 
data set shows that by the time the observation ended, 
the pulse profile had not recovered to its pre-burst shape. 

The first three bursts occurred near the middle “new” 
peak in the profile. Burst 5 occurred near the old tall 
peak of the profile (see Fig. 5). The phase of burst 6 
corresponds again to a new peak in the profile, this time 
located where the previous small peak used to be. The 
coincidences of several bursts with new, transient pro- 
file features that are present even when the actual burst 
data are removed are suggestive of lower-level transient 
emission from the same physical location, with the burst 
being the extreme of this emission’s luminosity distribu- 
tion. 


3.3. Burst Rotational Phases 

An important factor in understanding the origin of the 
bursts is the rotational phase at which they occur. The 
phases of the bursts are shown in Figure 7. For each 
burst observation we created at 31.25-ms time resolution 
lightcurve and folded it using our timing solution (see 
§ 3.5). We then phase-aligned these folded profiles by 
cross-correlating them with the long-term pulse profile 
template. Our phase- aligned folded profiles axe shown 
in Figure 7 (histograms in the last row). Superposed 
on each folded profile are two curves. The top curve is 


3.4. Pulsed Flux Analysis 
3.4.1. Short-Term Pulsed Flux Analysis 

Previous AXP bursts are often accompanied by short- 
term pulsed flux enhancements (e.g. Gavriil et al. 2006). 
To search for these, for each of the three observations 
containing bursts, we made two time series in count rate 
per PCU, one for the 2—4 keV band and the other for 
4—20 keV. We included only photons detected by PCUs 
that were on for the entire duration of the observation. 
The time resolution was 31.25 ms. We removed the 4 s 
centered on each burst from each time series. Then, we 
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Fig. 4. — Spectral evolution of burst 6. Top: Burst lightcurve binned with 4.34-s (half a pulse period) time resolution and extracted 
over the same energy band used for spectral fitting, 2.5-30 keV. Middle: Surface plot of the difference between the burst spectrum and 
a simple continuum model (blackbody) in terms of change in x 2 - Each vertical slice in the surface plot corresponds to a spectrum that 
was calculated over a two-pulse-period long window. This window was translated across the burst in half-a-puise-period long steps. For 
each spectrum, a best-fit blackbody model was subtracted. The color wedge on the right maps the colors to the resulting change in x> 
The contours indicate the regions where the spectrum differed by the model by 1 -a (blue contour), 2-er (green contour) and 3-<r (yellow 
contour), respectively. The feature at 14.2 keV is highly significant at the start of the burst, and remains detectable at the >1 -<j level for 
~130 s. Notice, that the significance of the 8.6 keV feature is intermittent, becoming comparable to that of the 14.2 keV feature at later 
times, and is the only feature detectable at the >l-<7 beyond ~130 s from the burst peak. Bottom: The insets are examples of individual 
two-pulse-period long spectra used to generate the surface plot. These spectra were extracted at 8.69, 52.1 and 139.0 s from the burst 
onset, respectively. The red curves indicate the corresponding best-fit blackbody model. 


broke each time series into segments of length ~500 s. 
For each segment, we calculated the pulsed flux using 
two different methods. 

First, we calculated the RMS pulsed flux using 


Frms 


n 


N 


2j2({a k 2 + b k ' 2 )- (<r ajl 2 

k=l 


+ <0), 


(4) 


where a* is the k even Fourier component defined as 
a k = ipiC0s(27rfcz/iV), cr ak is the uncertainty of 

a ky b k is the kfi 1 odd Fourier component defined as b k 
= ■^^^ 1 PiSin(27rfct/A/'), cr& fc is the uncertainty of b k , i 
refers to the phase bin, N is the total number of phase 
bins, pi is the count rate in the i*'* 1 phase bin of the 
pulse profile, and n is the maximum number of Fourier 
harmonics to be taken into account; here n — 5 for consis- 
tency with Dib et al. (2007) and Gonzalez et al. (2010). 

While least sensitive to noise compared to other pulsed 
flux measurement methods, the RMS method returns 


a pulsed flux number that is affected by pulse profile 
changes. To confirm our pulsed flux results, we also used 
an area-based estimator to calculate the pulsed flux: 


N 


^Area ~ ^ r ^ J Pi Pmin 


(5) 


i = 1 


where p m i n is the average count rate in the off-pulse phase 
bins of the profile, determined by cross-correlating with 
a high signal-to-noise template, and calculated in the 
Fourier domain after truncating the Fourier series to 5 
harmonics. The results are shown in Figure 8 for Fr ms 
(^A rea gives consistent results). Note the significant in- 
crease in the 4—20 keV pulsed flux in the 2006 June ob- 
servation following the cluster of bursts. This increase is 
not present in 2—4 keV. Also note the significant rise and 
subsequent decay of the pulsed flux following the large 
2007 February burst. The pulsed flux was sufficiently 
enhanced in these two observations that one can see in- 
dividual pulsations by eye in Figure 7 in the two panels 



Fig. 5. — Pulse profile evolution of 4U 0142+61. Left-most column: normalized average 2—10 keV pulse profiles in the 6 years prior to 
the entry into the active phase in chronological order (panels from Dib et al. 2007). Middle and right-most columns: normalized, averaged 
2—10 keV pulse profiles for each of the data groups marked by a letter in Figure 1, after the entry into the active phase. The 4 plots with 
bold labels correspond to the observations marked with “x”, “bl”, “b2’\ and “b3” in Figure 1. 


marked with a star, containing the raw burst lightcurves 
with 1-s time resolution. 

3.4.2. Long-Term Pulsed Flux Analysis 

For each of the 94 analyzed observations, we created a 
pulse profile (in units of count rate per PCU) using the 
same procedure as in Section 3.2. Then we calculated 
the pulsed flux for each observation using Equations 4 
and 5. Data from PCU 0 were omitted because the long- 
term trend in the pulsed counts is not the same as that 
in the remaining PCUs, presumably due to the loss of its 
propane layer. Data from PCU 1 were omitted after the 
loss of its propane layer as well, on MJD 54094. 

We extracted pulsed fluxes in the 2—4 and 4—20 keV 
bands using both the RMS and area pulsed flux method 
because of the numerous pulse profile changes around the 
times of the bursts. There is no evidence of a iong-term 
change in the pulsed flux associated with the bursts ex- 
cept, possibly, for a hint of an increase in the 2—4 keV 
band roughly midway between the second and third ob- 


servations containing bursts. We also note that the 
4—20 keV pulsed flux of two of the observations con- 
taining bursts are significantly larger than the long-term 
average. Removing the bursts from these observations 
does not change this result. 

We performed the same analysis for individual obser- 
vations in 2—10 keV for an extended period of time. This 
is shown in panels 2 and 3 of Figure 9. This plot is an 
update to Figure 10 of Dib et al. (2007). 

In order to verify that trends seen in panels 2 and 3 
of Figure 9 are not an artifact of the evolution of the 
response of the detector, we calculated the pulsed flux 
in erg s -1 cm -2 using a method that took the evolution 
of the response into account. The method is described 
in Dib et al. (2007) and takes into account spectral fits 
obtained from imaging data (in this case XMM data, 
see Gonzalez et al. 2010) to convert counts to energy for 
each combined set of observations. For this analysis, we 
used data from all PCUs; however, data from PCU 1 
were excluded after the loss its propane layer. Data from 



2002 2004 2006 



Fig. 6. — Top: evolution of the power in the fundamental Fourier component of the pulse profile of 4U 0142-1-61 in the 2—10 keV 
band. Bottom: Evolution of the power in the second harmonic. In both panels, the points with large horizontal error bars are from Dib 
et al. (2007). The remaining points are obtained from groups of observations after the entry into the active phase. The three bold points 
correspond to the observations with bursts. The three dashed lines correspond to burst epochs. 


PCU 0 were included because the response matrices used 
took into account the loss of its propane layer. 

The results are shown in panel 1 of Figure 9. The 
first 6 points, corresponding to RXTE Cycles 6—10, are 
from Dib et al. (2007). The second to last point is ob- 
tained by combining all observations that occurred be- 
tween bursts 1 and 2, but omitting observations contain- 
ing bursts. The observations that we included took place 
during the exponential recovery of the possible glitch (see 
Section 3.5). The last point in panel 1 was obtained by 
combining the observations that occurred after burst 2, 
but again omitting those containing bursts. The obser- 
vations we included took place after the end of the expo- 
nential recovery of the possible glitch (see Section 3.5). 

In the first of the two data points in the active phase, 
the pulsed flux in erg s _1 cm -2 is 18+8% larger than 
in the pre-active phase. This is consistent with the in- 
crease reported in Gonzalez et al. (2010) in the same 
energy range. A hint of this increase can be seen in 
panel 2 although it appears less significant. This discrep- 
ancy can be accounted for by the fact that the spectrum 
changed: Gonzalez et al. (2010) reported a temporary in- 
crease in the spectral hardness in an XMM observation 
of 4U 0142+61 immediately following the bursts. 

3.4.3. Combined Pulse Morphology and Pulsed Flux Analysis 


In Section 3.2, we calculated the Fourier components 
of the aligned average pulse profiles. This gave us the 
relative amplitude of the pulse profile harmonics in each 
group of observations marked with a letter in Figure 1. 
In Section 3.4.2, we calculated the pulsed flux for every 
observation. Here, we compute a weighted average of 
the pulsed flux for each group of observations using the 
flux points calculated in Section 3.4.2. We then recon- 
struct the profiles for each of the groups from the first 
five Fourier components (not including the DC), scale 
them by the average RMS pulsed flux for that group, 
and add the necessary offset for the lowest point on each 
curve to be zero. This means that the resulting scaled 
profiles return the correct pulsed flux. The advantage of 
this analysis is that we can now trace the evolution of 
each of the peaks independently. 

The results are presented in Figure 10. In panel (a), we 
show the average scaled pulse profile per observation for 
observations in the year preceding the active phase. The 
profile is double-peaked. In panel (b), we show in bold 
the scaled profile for the observation containing burst 1 . 
We also show the scaled profiles in each of the groups 
of observations that followed the burst. The profiles are 
triple-peaked. In panel (c), we show in bold the scaled 
profile for the observation containing bursts 2, 3, 4, and 5. 
We also show the scaled profiles in each of the groups of 
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Fig. 7. — Phase analysis of the 6 detected bursts. Bach column corresponds to an observation in which bursts were detected. Tbp row: 
200-s segments of the time series containing the bursts in the 2—20 keV band. The time resolution is 1 s. Middle row: 4-s segments of 
the time series containing the bursts. The time resolution 31.25 ms. Bottom row: aligned folds of the burst observations shown below 
the scaled long-term average profile. The folded counts are presented as histograms. Superposed on each histogram are two curves. The 
top curve is made of the 5 Fourier harmonics that best fit the histogram. The bottom curve is made of the best-fit 5 harmonics after the 
removal of the 4 seconds centered on each burst. The arrows indicate the phases of each burst. 
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Fig. 8. — RMS pulsed flux within the observations containing bursts. Each column corresponds to one observation. In each column we 
show, descending vertically, the 1-s resolution lightcurve with the bursts indicated, the 2-4 keV RMS pulsed flux, and the 4-20 keV RMS 
pulsed flux. The dotted line in each of the pulsed flux plots shows the average of the pulsed fluxes obtained after segmenting and analyzing 
the time series of the observation immediately prior to the one shown. 
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Fig. 9. — Long-term pulsed flux time series in 2—10 keV 
for 4U 0142+61. Panel 1: area pulsed flux in erg s -1 cm" 2 
for combined observations. The solid line marks the 29+8% in- 
crease reported in Dib et al. (2007). Panel 2: area pulsed flux in 
counts /s/PCU for individual observations. Panel 3: RMS pulsed 
flux in counts/s/PCU for individual observations. All panels: the 
dashed lines mark the burst epochs. The points marked with stars 
are the pulsed fluxes of the observations containing bursts. 

observations that followed the bursts. Note how the in- 
crease in the pulsed flux in the observation containing the 
bursts is not only a consequence of the appearance of the 
new peak, but a result of the increase in size of all three 
peaks. Also note the evolution of the pulse profile back 
to being double-peaked. In panel (d), we show in bold 
the scaled profile for the observation containing burst 6. 
We also show the scaled profiles in each of the groups 
of observations that followed the burst. Note again how 
the pulsed flux increase is due to both peaks increasing 
in size. 


3.5. Timing Analysis 

Many AXP outbursts and active phases are accompa- 
nied by interesting timing anomalies (e.g. Kaspi et al. 
2003; Woods et al. 2004; Israel et al. 2007a; Dib 
et al. 2009). Next we consider the timing behavior of 
4U 0142+61, which, as we show, also exhibits interesting 
evolution at the start of the active phase. 

For all our RXTE observations of 4U 0142+61 we cre- 
ated 31.25-ms time resolution lightcurves, including only 
the events in the energy range 2.5—9 keV to maximize 
the signal-to-noise ratio of the pulse. Each binned time 
series was epoch-folded using an ephemeris determined 
iteratively by maintaining phase coherence; see below. 



Phase (Cycles) 


Fig. 10. — Average pulse profile per observation in the 2—10 keV 
band in several groups of observations reconstructed from the first 
five Fourier components, and scaled to return the appropriate 
pulsed flux, (a) average scaled pulse profile for observations in 
RXTE Cycle 10, the last RXTE Cycle before the active phase. 
(b) Bold curve: scaled pulse profile of the observation containing 
burst 1. Thin curves: average scaled pulse profile for each of the 
groups of observations that followed burst 1 (groups G and H in 
Figure 1). (c) Bold curve: scaled pulse profile of the observe 

tion containing bursts 2, 3, 4, and 5. Thin curves: average scaled 
pulse profile for each of the groups of observations that followed 
the bursts (groups I, J, K, and L in Figure 1). (d) Bold curve: 
scaled pulse profile of the observation containing burst 6. Thin 
curves: average scaled pulse profile for each of the groups of ob- 
servations that followed burst 6 (groups M, and N from Figure 1). 
The vertical dotted lines indicate the phases of each burst. 


Resulting pulse profiles, with 64 phase bins, were cross- 
correlated in the Fourier domain with a high signal-to- 
noise template created by adding phase-aligned profiles 
from all observations. The cross-correlation returned an 
average pulse time of arrival (TOA) for each observation 
corresponding to a fixed pulse phase. The pulse phase 0 
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Fig. 11. — Top: Timing residuals for all three sets of TOAs. Bold circles: timing residuals obtained using the first set of TO As obtained 
by cross correlation in the Fourier domain. Empty circles: residuals obtained using the second set of TOAs obtained by aligning the 
off-pulse regions of the pulse profiles. Triangles: residuals obtained using the third set of TOAs obtained by aligning the tallest peak of the 
pulse profiles. The linear ephemeris used to produce all three sets of residuals is that which best fit the first set of TOAs in the pre-active 
phase region. Bottom: Timing residuals for the same three sets of TOAs. The linear ephemeris used is that which best fit the first set of 
TOAs exclusively in regions indicated with the horizontal arrows with an arbitrary phase jump in the gap between the arrows. The phase 
residuals are relative to the linear ephemeris (ignoring the phase jump) and vertically offset such that the average phase residual is zero. 
Both panels: the dotted line indicates the epoch of the first observation of the active phase. The dashed lines indicate the burst epochs. 


at any time t can be expressed as a Taylor expansion, 

<!>(t) = <j>o(to) + Vo(t - to) + 2^°^ — ^°) 2 

+ ^ va(t — to) 3 + • • «i (6) 

where v = 1/P is the pulse frequency, 0 = dv/dt , etc., 
and subscript 4 0” denotes a parameter evaluated at the 
reference epoch t = to- To obtain ephemerides, the TOAs 
were fitted to the above polynomial using the pulsar tim- 
ing software package TEMPO 7 . 

As explained above, our first set of TOAs was obtained 
by aligning the folded observations with a template pro- 
file using a cross-correlation procedure. In order to de- 
termine to what extent this set of TOAs was affected by 
the pulse profile changes that took place in the active 
phase, we generated two additional sets of TOAs (sets 2 
and 3). 

For set 2, we made the assumption that the location 
of the trough of the pulsations, determined by finding 

7 See http: //www. atnf.csiro.au/research/pulsar/tempo. 


the minimum in the pulse profile after smoothing, is not 
affected by the pulse profile changes. We then generated 
TOAs by aligning the pulsation trough of the folded ob- 
servations with that of the template and extracting phase 
differences. The last set of TOAs (Set 3) was obtained 
after aligning the tallest peak of each smoothed pulse 
profile with that of the long-term template. All three 
sets of TOAs are provided as a machine readable table 
in the electronic edition. A portion of the table starting 
from the active phase is provided in Table 3. 

Timing residuals for all three sets of TOAs axe shown 
in Figure 11 using two different ephemerides. The linear 
ephemeris used in the top panel is the ephemeris that 
best fit the first set of TOAs prior to the active phase. 
The linear ephemeris used in the bottom panel is that 
obtained by fitting the first set of TOAs in the regions 
indicated by the horizontal arrows, with a phase jump 
in between. The same long-term trend is present in the 
first two sets of timing residuals, with more scatter in the 
second set. The difference in phase between each bold 
circle and the corresponding empty circle represents our 
uncertainty in determining a fiducial point on the pulsar. 
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TOA 

Numb.® 

Obs. ID b 

Exp. c 

(ks) 

Active 

PCUs d 

Set 1 

TOA e 

(MJD) 

°TOA f 

(s) 

Set 2 
TOA® 

(MJD) 

tf-TOA 5 

(•) 

Set 3 
TOA® 

(MJD) 

CTOA f 

(s) 

PT 

Analysis?® 

33 

92006-05-02-00 

4.9 

3 

53817.1615364 

0.08 

53817.1616238 

0.20 

53817.1615403 

0.20 

Yes 

34 

92006-05-03-00 

4.5 

3 

53831.3358686 

0.09 

53831.3358558 

0.21 

53831.3358711 

0.21 

Yes 

58 

92006-05-09-01 

1.0 

3 

53911.0530210 

0.17 

53911.0531137 

0.43 

53911.0530298 

0.43 

Yes 

78 

92006-05-25-00 

5.3 

2 

54138.4487229 

0.11 

54138.4487002 

0.29 

54138.4487277 

0.29 

No 


“ TO A number. TOA number 33 is the first observation in the active phase. TOAs number 34, 58, and 78 represent the first, second 
and third observation containing bursts. 
b RXTE observation identification number. 
c The effective number of PCUs on during the observation. 
d Total exposure. 

e Barycentered pulse TOA. Set corresponds to the TOAs obtained by cross- correlating the folded pulse profiles with a high signal- 
to-noise template. Sets 2 and 3 are similar to Set 1, however, we constrained our cross-correlation to the trough and peak of the 
pulsations, respectively. See § 3.5 for details. 
f Uncertainty on TOA for Set 1, 2, and 3. 

g Whether the TOA was included in the partial timing analysis. Note: All TOAs were included in the long-term analysis. 


In the third set of residuals, the outlier points represent 
the observations where the largest peak was no longer 
the right-most peak. Apart from the outliers, the same 
trend seen in the other two sets of residuals is present in 
the third set. 

From here on, we assume that the presence of the same 
trend in our three sets of residuals is an indication that 
the TOAs in the first set were not significantly affected 
by the pulse profile changes. We therefore have used the 
first set of TOAs (one for each of the 94 observations) in 
the remainder of this Section. 

Using the above assumption, the results of the tim- 
ing analysis are shown in Figure 12. The results of a 
segmented timing analysis are in panels 1 and 2. The 
results of a long-term timing analysis are in panels 3, 4, 
5, and 6 of the same Figure. 

For the segmented timing analysis, we used the 70 
TOAs that had the smallest uncertainties (indicated in 
the online version of Table 3) and omitted the TOAs of 
the 3 burst observations as well as the TOA of the first 
observation in the active phase. We then divided the 
data into segments of similar pulse profiles. For each data 
segment we found a linear ephemeris using TEMPO 8 . We 
plotted the resulting ephemerides in panel 1 of Figure 12 
with the uncertainties returned by TEMPO. The timing 
residuals are shown in panel 2. Notice how the slope be- 
tween bursts 1 and 2 is more negative than the long-term 
average. 

We then fit a linear trend through all the pre-burst 
observations, and another linear ephemeris through the 
post-burst observations after the pulse profile had started 
to return to the double-peaked shape. We did not include 
the TOA corresponding to the last observation contain- 
ing bursts. The results are plotted in panel 3 of Figure 12 
with uncertainties, with the residuals in panel 4. Notice 
the difference in the slope in the two regions. In partic- 
ular, extrapolating the two ephemerides to a point be- 

8 Fbr the segmented timing analysis we ran TEMPO in mode 1. 


tween their times of validity makes it seem like a sudden 
spin down, i.e. an “anti-glitch” occurred. 

Finally, we included all the TOAs and did one global 
fit. In order to provide a good fit to the TOAs at the 
onset of the active phase, we had to assume that a glitch 
occurred on MJD 53809, with the glitch model consisting 
of a permanent change in z/ and z> and a frequency change 
Ai/d that recovered exponentially on a timescale r^, i.e., 

V - U 0 {t) + Av + Av i e- { - t - t ’ )/Ti + AO (t - t g ), (7) 

where i/q (t) is the frequency evolution pre-glitch given 
by vo(t) = v(to) + i/(t — £o)> Az/ is a instantaneous un- 
recovered frequency jump, Ai/^ is the frequency increase 
that decays exponentially on a time scale r^, t g is the 
glitch epoch determined by setting the phase jump to 
zero, and A z> is the post-glitch change in the long-term 
frequency derivative. The values of the fit parameters 
are listed in Table 4. 

From Table 4, we can see that the total sudden change 
in frequency which happened at the onset of the active 
phase (A u tot = Az/ 4- A Vd) is positive, and A v& decayed 
exponentially. At the end of the decay, the net effect 
of the recovered glitch was a negative Av. That is, the 
data suggest that after v& had decayed, there remained 
net spin -down relative to the undisturbed ephemeris from 
before the active phase. 

The fit parameter Az/<*, from which we conclude a sud- 
den spin-up, is controlled primarily by the first few TOAs 
in the active phase. It is therefore possible that it is af- 
fected by pulse profile changes. However, as shown in 
Figure 11, the TOAs in the active phase were not signif- 
icantly contaminated by pulse profile changes. Also, the 
segmented analysis shown in Figure 12 clearly indicates 
an initial spin-up. 

In contrast to Az/<j, the fit parameter A z/, from which 
we conclude a net spin down , is primarily controlled by 
the TOAs outside of the active phase, which were cer- 
tainly not affected by pulse profile changes (Figure 12, 
panel 3). 
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Fig. 12. — Panel 1: Frequency vs. time plot of local linear 
ephemerides with uncertainties. The linear trend corresponding to 
the local ephemeris directly before the active phase is subtracted 
from all local ephemerides. Panel 2: Timing residuals after sub- 
tracting the ephemerides shown in panel 1. Panel 3: Frequency vs. 
time plot (with uncertainties) of the best-fit linear ephemerides in 
the pre-active phase region and in the post burst 2 region. The 
same trend as in panel 1 was subtracted. Panel 4: Timing resid- 
uals after subtracting the ephemerides shown in panel 3. Panel 5: 
Frequency vs. time plot (with uncertainties) of the best-fit long- 
term ephemeris which includes a glitch at MJD 53809 followed by 
a fast exponential recovery. The same trend as in panel 1 was sub- 
tracted. Panel 6: Timing residuals after subtracting the ephemeris 
shown in panel 5. All panels; the dotted line indicates the epoch of 
the first observation of the active phase. The dashed lines indicate 
the burst epochs. 


Parameter 

Value 

MJD start 

53481.268 

MJD end 

54235.319 

TOAs 

93 b 

is (Hz) 

0.1150920955(12) 

is (10 -14 Hz s _1 ) 

-2.661(9) 

Epoch (MJD) 

53809.185840 

Glitch Epoch (MJD) 

53809.185840 

Ai/ (Hz) 

— 1.27(17)xl0 -8 

A is ( Hz s -1 ) 

— 3.1(1. 2)xl0 -16 

A v d (Hz) 

2.0(4) xl0~ 7 

t d (days) 

17.0(1.7) 

RMS residual (phase) 

0.0168 


a Numbers in parentheses are TEMPO-reported lcr uncertainties. 
b A single TO A was omitted due to the very poor ‘signal to noise 
ratio in the corresponding observation. 

We therefore conclude that the pulsar likely suffered a 
spin-up glitch near or at the start of the radiatively active 
phase, but that the glitch ‘over-recovered’ such that long- 
term, its net effect is a spin -down of the pulsar. 

4. DISCUSSION 

In this paper, we have described the timing, pulse pro- 
file, and pulsed flux behavior of 4U 0142+61, during its 
2006-2007 active phase, the first such episode yet stud- 
ied from this source. Specifically we have shown that in 
addition to a sudden departure from a slow, systematic 
evolution of the source’s pulsed flux and pulse profile, 
this AXP also suffered a significant timing event that is 
best described as a sudden spin-up glitch, followed by 
a large decay of the frequency jump such that the net 
effect was a slow -down with Av/is ^ -8 x 10 -8 . In- 
terestingly, the pulsed X-ray flux showed no significant 
change apart from immediately following bursts, in con- 
trast to other AXP radiative outbursts (e.g. Kaspi et al. 
2003; Tam et al. 2008; Ibrahim et al. 2004; Muno et al. 
2007). Further, we have reported on six SGR-like bursts 
from the source that occurred during this period, one of 
which was notable for its unusual spectrum, which was 
poorly fit by an continuum model. 

4.1. The SGR-Like Bursts from 4U 0142+61 

Five of the six bursts reported here for 4U 0142+61 had 
fast-rise-slow-decay profiles, with tails much longer than 
the rise times. These morphologies are similar to the 
class of AXP bursts labeled as “Type B” by Woods et al. 
(2005). Type B AXP bursts have also been character- 
ized by arriving preferentially at pulse maximum, and are 
seen more often in AXPs compared with SGRs (Woods 
et al. 2005; Gavriil et al. 2006). They were suggested. as 
being due to a sudden rearrangement of magnetic field 
fines anchored in the crust following a crustal fracture 
(Thompson & Duncan 1995), as opposed to reconnection 
events in the upper magnetosphere (e.g. Lyutikov 2002), 
argued as more likely for the shorter, symmetric “Type 
A” bursts that show no preference for pulse maximum. 
However, it is interesting that in the 4U 0142+61 bursts, 
in spite of their B-type morphology, no clear preference 
for arrival at or near pulse maximum was seen (see Table 
1). This may blur somewhat the distinction between the 
putative types, although the morphological distinction 
remains clear. 
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the minority (roughly one dozen out of 80) of the bursts 
seen were of this form, the majority being of Type A, 
similar to those classically seen in SGRs, and suggested 
by Woods et al. (2005) to be magnet ospheric. A possible 
hint regarding the origin of the different types may lie 
in that IE 2259-1-586 was mid-outburst when its bursts 
were observed, whereas the other sources* bursts all oc- 
curred in the days/months following the commencement 
of an active period, presumably as the pulsar recovered 
from a major event, rather than mid-event. 

4.1.1. Spectral Features 

The complicated spectrum of burst 6 is puzzling. The 
most broad feature at ~14 keV is particularly interesting. 
Emission features at similar energies were observed from 
two out of the three bursts from IE 1048.1-5937 (Gavriil 
et al. 2002; Gavriil et al. 2006) and in one out of the four 
bursts from XTE J1810-197(Woods et al. 2005). We 
have reanalyzed all these burst spectra consistently as 
for 4U 0142-1-61. In Figure 13 w r e plot the spectra of all 
AXP bursts with likely emission lines in their spectra. 
Notice that all have broad features that occur between 
13 and 14 keV. There is also evidence for a narrow feature 
at ~8 keV in the 4U 0142+61 burst; there may be hints 
of additional features in some of the other burst spectra 
as well. 

If the features are interpreted as of a cyclotron origin, 
it is puzzling that three sources with spin-down-inferred 
magnetic fields that span a range of a factor of three 
produce features at similar energies. The similarity in 
energy suggests a mechanism that is not very sensitive 
to the exact value of the magnetic field. Moreover, it is 
unclear why such features are seen only in some bursts. 
Cyclotron features have been observed in the persistent 
spectra of accretion-powered X-ray pulsars (Coburn et al. 
2002). Unlike the features observed here, these features 
can be well fit by a Gaussian absorption model (see 
§ 3.1.2). In the spectra of some accretion-powered X- 
ray pulsars there is also evidence for higher harmonics. 
The features observed in the 4U 0142+61 burst were not 
consistent with being harmonically related, even after 
taking into account the relativistic effects on the spacing 
(see Harding k Daughtery 1991; Araya & Harding 1999). 

Thus far, all AXP burst spectral features have been 
observed with RXTE. Their absence from many other 
#XT#-observed bursts, together with the absence of any 
response feature in the PCA near 13-14 keV, is support- 
ing evidence for their veracity. Nevertheless, it is impor- 
tant to confirm these features with other instruments. 
Observations of AXPs during active phases using the 
Swift XRT or the large-area LAXPC aboard the planned 
Astrosat mission would be particularly useful in this re- 
gard. 

4.2. The Net Spin-Down Event in 4U 0142+61 

The timing glitch reported on in §3.5 had recovery 
fraction, defined as Q = Avd/{Av d + Ap), equal to 
1.07 + 0.02. Q > 1 implies that the net frequency change 
after the transients have decayed is negative (see panel 3 
of Figure 12 and Table 4). Indeed, the net Av/v ob- 
served is -1.1 ± 0.1 x 10“ 7 . Such a Q > 1 has not 



Fig. 13. — Burst spectra of all observed AXP bursts with signifi- 
cant emission lines, as seen by RXTE. Panel 1A: Burst spectrum of 
4U 0142+61 burst 6. The solid line indicates the continuum (black- 
body) component of the best-fit model folded through RXTE J a in- 
strumental response. Panel IB: Residuals, expressed in terms of 
the change in x 2 > after subtracting only the continuum component 
of the best-fit model. Panels 2 A and 2B: Same but for burst 4 of 
XTE Ji 810-197 (see Woods et al. 2005). Panels 3 A and 3B: Same, 
but for burst 3 of IE 1048.1-5937 (see Gavriil et al. 2006). Panels 
4A and 4B: Same, but for burst 1 of IE 1048.1-5937 (see Gavriil 
et al. 2002). 

previously been reported in any AXP but has been seen 
following magnetar-like radiative behavior in one high-H 
rotation-powered pulsar (Livingstone et al. 2010). Note 
that Morii et al. (2005) reported a candidate glitch in 
1999 from this source. If such a glitch occurred, based 
on the ephemerides reported in Dib et al. (2007), the 
fractional frequency change, after any recovery, would 
have been +1.9(2) xlO -7 — 8.6(2) xl0“ 7 depending on the 
glitch epoch, but manifestly positive. 

In the standard model for glitches in rotation-powered 
pulsars, the neutron-star crust contains superfluid neu- 
trons that rotate faster than the bulk of the surrounding 
matter (see, e.g., Anderson & Itoh 1975; Anderson et al. 
1982; Alpar et al. 1984; Alpar & Pines 1993). This is ar- 
gued to be a result of the fact that the external magnetic 
torque acts on the crust and coupled core components 
only, with the uncoupled superfluid unaffected. The su- 
perfluid's angular momentum resides in quantized vortex 
lines whose density is proportional to angular velocity. 
The vortex lines are suggested to be pinned to crustal nu- 
clei, and suffer strong forces due to the angular velocity 
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superfluid to crust. In a magnetar, unpinning may occur 
due to the strong internal magnetic field as it deforms 
the crust plastically or cracks it violently (Thompson & 
Duncan 1996). 

In the 4U 01424-61 event we describe in §3.5, if the 
standard glitch model is roughly correct and applies here, 
then some regions of the stellar superfluid were origi- 
nally spinning slower than the crust. Then the tran- 
sient increase in frequency would be a result of transfer 
of angular momentum first from a more rapidly rotat- 
ing region, with a subsequent angular momentum drain 
from the crust to a more slowly rotating region. As 
argued by Thompson et al. (2000), regions of slower- 
rotating superfluid can occur in magnetars, because the 
superfluid vortex motions are governed not be spin-down- 
related forces, but by aavection across the stellar surface 
by the deforming crust. Those authors show that the 
number of vortex lines per unit surface area of crust 
can increase or decrease depending on the crustal mo- 
tion relative to the stellar rotation axis. They invoke 
this possibility to explain a possible “anti-glitch” seen in 
SGR 19004-14. However, that event was orders of magni- 
tude larger than what we have observed in 4U 01424-61, 
with ~ 10 -4 (Thompson et al. 2000), and also was 
likely associated with the extremely energetic flare seen 
on 1998 August 27 (Mazets et al. 1999). 

Independent evidence for slow crustal deformations in 
4U 01424-61 that could result in the regions of slower- 
rotating superfluid required to explain the spin-down 
may come from the long-term flux and pulse profile evo- 
lution previously reported for this source by Dib et al. 
(2007) and shown again here in Figures 5 and 6. Al- 
though such variations could also be magnetospheric in 
nature, possibly due to twisting of the magnetic field 
(e..g. Thompson et al. 2002; Beloborodov &; Thompson 
2007), problems with this interpretation exist, as dis- 
cussed by Dib et al. (2007). For example, if magneto- 
spheric, it is puzzling that most if not all of the changes 
in pulse profile are seen below 4 keV, with none above 
6 keV. On the other hand, slow evolution of the sur- 
face emission (e.g. Ozel & Giiver 2007), hence structure, 
is perhaps more consistent with low-energy evolution. It 
remains to be seen if detailed modeling can reproduce the 
sort of slow pulse profile changes we have observed (Fig. 
5), quite apart from the sudden appearance of new peaks, 
hence greater harmonic structure, near bursts. We note 
that the latter has also recently been seen during periods 
of activity in IE 1048.1-5937 (Tam et al. 2008). 

Previously Kaspi et al. (2000) noted that the existence 
of AXP glitches having properties comparable to those 
seen rotation-powered pulsars was not in itself evidence 
for AXPs being magnetars. This is because the stan- 
dard glitch model applies regardless of the mechanism 
by which the crust slows down; i.e., glitches can occur, 
in principle, in an accretion-powered pulsar. For exam- 
ple, a glitch was a plausible explanation for the fractional 
change of A v/v ~ 4 x 10 -5 seen in the accreting neu- 
tron star system KS 19474-300 (Galloway et al. 2004). 
However, although glitches can occur in the case of ac- 
creting systems, since there too crust /superfluid angular 
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down in 4U 01424-61 may thus add to the already large 
amount of evidence against any accretion-powered origin 
for 4U 0142+61. 

4.3. The Active Phase of 4U 0142+61 and Other AXP 

Outbursts 

In many ways, the 2006-2007 active phase of 
4U 0142+61 is similar to other phases of activity seen 
in AXPs: it was punctuated by short bursts, pulse pro- 
file changes were seen, and it was accompanied by a sig- 
nificant rotational anomaly. However the 4U 0142+61 
activity is unique in one interesting way: the pulsar did 
not suffer a large, long-lived pulse flux increase at the 
beginning of the phase. Indeed its pulsed flux (Fig. 9), 
apart from very near bursts, has remained relatively sta- 
ble. We note that observations with focusing X-ray tele- 
scopes may reveal some phase-averaged flux variations 
(see Gonzalez et al. 2010); indeed pulsed fraction has 
been shown to be inversely correlated with total flux in 
IE 1048.1-5937 (Tiengo et al. 2005; Tam et al. 2008), 
rendering pulsed flux variations smaller relative to to- 
tal variations. Still, for 4U 0142+61, such a correlation 
would have to conspire to render the pulsed flux steady. 
This seems unlikely, however only focusing telescope ob- 
servations can rule this out. 

Other AXP radiative outbursts have been suggested 
as arising from large magnetospheric twists, with associ- 
ated magnetospheric currents returning to the stellar sur- 
face and heating it, resulting in increased X-ray emission 
from the source (e.g. Thompson et al. 2002; Beloborodov 
&: Thompson 2007; Beloborodov 2009). Such twists are 
thought to represent a release of magnetic energy and 
helicity from the internally wound-up magnetic field. 

On the assumption that the pulsed flux is a reasonable 
proxy for the total flux for 4U 0142+61, here, we find no 
clear evidence for significantly increased X-ray emission 
on time scales longer than a few minutes. Thus any large- 
scale magnetospheric twist scenario is problematic. As 
originally showed by Thompson et al. (2002), the X-ray 
luminosity due to returning currents reheating the sur- 
face in a significantly twisted magnetosphere in general 
are comparable to that produced from internal processes; 
this is clearly not observed in the 2006-7 active period of 
4U 0142+61. 

Instead, for this source, long-term evolution of the 
crust, driven presumably internally by field decay, and re- 
sulting in multiple unstable configurations though with- 
out any large scale magnetospheric twists, could result 
in sudden cracking and local rearrangements. This could 
be accompanied by bursts and profile changes, as well as 
with vortex line shifting. Why such surface motion does 
not produce significant magnetospheric twists is, how- 
ever, puzzling, given that field lines are thought to be 
anchored in the crust; perhaps only for large motions do 
field lines become sufficiently twisted for enhanced X-ray 
emission to be produced. 

Dib et al. (2008) showed that other AXP glitches have 
been unaccompanied by radiative changes. For exam- 
ple, AXP IE 1841-045 has glitched multiple times yet its 
pulsed flux remains very steady. Such radiatively “silent” 
glitches may be a result of internal activity that does 
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in fact that thus far, the data are consistent with all 
AXP radiative outbursts being accompanied by timing 
anomalies. This may also be true of SGRs, though in 
those cases, timing anomalies are harder to establish be- 
cause of the difficulty in achieving phase-coherent timing 
in quiescence. 

5. SUMMARY AND CONCLUSIONS 

We have reported on anomalous X-ray pulsar 
4U 0142—61 entering an active phase which was preceded 
with by a long-term increase in pulsed flux. The active 
phase, which commenced in 2006 March, consisted of a 
timing anomaly that can be described as a net “anti- 
glitch,” that is, a net spin-down following an initial spin- 
up that decayed on a time scale of 17 days. Following 
the glitch, we detected six bursts from the pulsar, the 
first ever observed from this source. Despite 10 yr of 
RXTE monitoring, the bursts all occurred in the nar- 
row time span between 2006 April 6 and 2007 Febru- 
ary 7. The sixth and largest burst had a unusual spec- 
trum that cannot be described by any simple continuum 
model. Rather, it can be well fit if a broad spectral 
feature at ~14 keV is included, as has been reported in 
other AXP bursts, as well as a narrow feature at ~8 keV. 
The pulse profile of the source changed from double- to 
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pulse profile is now relaxing to its pre-active phase mor- 
phology. Most aspects of 4U 0142+61’s emission changed 
during the active phase, with the notable exception of the 
pulsed flux (except near bursts). This argues against this 
event being associated with a sudden magnetospheric 
twist, as has been invoked for other AXP activity, and 
is suggestive of crustal evolution driven internally by the 
large magnetic field, though without significant magne- 
tospheric twisting. We suggest that other, radiatively 
silent AXP glitches have a similar origin, whereas radia- 
tively loud AXP timing events occur when the crustal 
motions cannot avoid significant twisting of the magnetic 
field lines. 
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